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Abstract. We investigate the interacting dark energy models by using the diagnostics of 
statefinder hierarchy and growth rate of structure. We wish to explore the deviations from 
ACDM and to differentiate possible degeneracies in the interacting dark energy models with 
the geometrical and structure growth diagnostics. We consider two interacting forms for 
the models, i.e., Qi = [3Hpc and Q 2 = with /3 being the dimensionless coupling 

parameter. Our focus is the lACDM model that is a one-parameter extension to ACDM by 
considering a direct coupling between the vacuum energy (A) and cold dark matter (CDM), 
with the only additional parameter /3. But we begin with a more general case by considering 
the IrcCDM model in which dark energy has a constant w (equation-of-state parameter). 
Eor calculating the growth rate of structure, we employ the “parametrized post-Eriedmann” 
theoretical framework for interacting dark energy to numerically obtain the e(z) values for 
the models. We show that in both geometrical and structural diagnostics the impact of w is 
much stronger than that of /3 in the IrcCDM model. We thus wish to have a closer look at 
the I ACDM model by combining the geometrical and structural diagnostics. We find that 
the evolutionary trajectories in the plane exhibit distinctive features and the depar¬ 

tures from ACDM could be well evaluated, theoretically, indicating that the composite null 
diagnostic is a promising tool for investigating the interacting dark energy models. 

We also compare our results with the observed uncertainties on diagnostic parameters. We 
find that current observations still do not have sufficient precisions to completely distinguish 
lACDM models from the ACDM model. Anyway, our work points out what precisions of 
measurements should be achieved to distinguish the lACDM models from the ACDM model. 
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1 Introduction 

The cosmological observations [1-4] have shown that our universe is undergoing accelerated 
expansion. This cosmic acceleration is attributed to an unknown component, called dark 
energy (DE), with negative pressure, or a modification to general relativity (GR) on the 
cosmological scale. For the possible existence of DE, a large number of DE models are 
constructed in theory, such as the ACDM model consisting of the cosmological constant 
(A) and the cold dark matter (CDM), the scalar-held models [5-15], the holographic type 
models [16-22], and so on. On orientation of modihed gravity, there are many attempts, for 
instance, the /(R) class theories [23, 24], the Dvali-Gabadadze-Porrati (DGP) braneworld 
model [25], and the like. Among a lot of theories and models, the ACDM model is the simplest 
and provides a good ht to current observational data, though the other models cannot be 
completely excluded yet. Nevertheless, the ACDM model is not free of any troubles; it 
always suffers from the severe theoretical challenges, i.e., the hne tuning and coincidence 
problems [26-32]. 

Currently, in the empirical studies of cosmic acceleration, we are actually facing two 
physically profound questions; Is the cosmic acceleration caused by a breakdown of GR on 
cosmological scales or by a DE component that produces repulsive gravity within GR? If 
the cause of acceleration is indeed a DE component, is it a cosmological constant A or a dy¬ 
namical field? To differentiate DE and modihed gravity (MG), the only usage of geometrical 
diagnostics is not enough, and the utilization of probes of structure growth is a must. In 
this paper, we do not consider the possibility of MG, but only focus on the studies of DE. 
But, actually, even if only the possibility of DE is considered, the combination of geometrical 
and structure growth diagnostics is proven to be very powerful and helpful in differentiating 
alternative models from AGDM, as will discussed in this paper. 

Undoubtedly, any dehnitive deviation from ACDM, if diagnosed out and conhrmed, 
would be a major breakthrough in cosmology and fundamental physics. One possible devia¬ 
tion comes from the models in which the cosmological constant A is replaced by a dynamical 
dark energy. There exists, however, another possibility that DE (in the simplest case, the 
vacuum energy density or A) directly couples to CDM. The model of vacuum energy interact¬ 
ing, in some physically profound way, with CDM is called the lACDM model in this paper. 
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We will explore the deviation of this model from ACDM with the combination of geometrical 
and structure growth diagnostics in this work. 

To discriminate various models, the geometrical diagnostics are introduced, such as 
statefinder {r, s} [33, 34], Om, and Om3 [35-37]. Some works [38-47] revealed that the 
statefinder diagnostic can effectively discriminate various DE models and break/differentiate 
the degeneracies between different parameter values of DE model, except for a few mod¬ 
els (for example, the new agegraphic DE model [45, 47]). Furthermore, Arabsalmani and 
Sahni [48] introduced the statefinder hierarchy and the growth rate of linear perturbations as 
“null diagnostics” for the ACDM model. The statefinder hierarchy [48] is also a geometrical 
diagnostic and is model-independent, which contains high derivatives of a(t), with a{t) the 
scale factor of the universe. The growth rate of structure [50, 51] was presented previously 
as a scale-independent consistency check between the expansion history and the structure 
growth. It can act as a cosmic growth history diagnostic, or be combined with the statefinder 
hierarchy to serve on a composite diagnostic. In the previous work [49], we applied these 
two diagnostics to discriminate four holographic DE models. These diagnostics were also 
considered in [52-54]. 

In this paper, we will use the statefinder hierarchy and the growth rate of structure 
to diagnose the deviation from the ACDM regarding the coupling between vacuum (A) and 
CDM in the lACDM model. Though we pay more attention to the lACDM model, we wish 
to begin with the more general model of interacting DE in which the equation of state (EoS) 
parameter w is a constant and may not be exactly —1, usually referred to as the IrcCDM 
model. We will first apply the geometrical and structural diagnostics in the ItcCDM model, 
exploring and comparing the impacts of EoS w and coupling /3. 

Consider a spatially flat Friedmann-Robertson-Walker (FRW) universe consisting of DE 
(de), CDM (c), baryons (b), and radiation (r). The energy conservation equations of cosmic 
components are expressed as 


Pde “1“ Q, 

(1.1) 

Pc “1“ Pc — Q J 

(1.2) 

Ph + = 0, 

(1.3) 

Pr-\-4:Hpr = 0, 

(1.4) 


where the dot denotes a derivative with respect to time t, H = d/a is the Hubble parameter. 
Pi is the energy density of each component, for i = de, c, b, and r, respectively, and Q is the 
energy transfer rate between DE and CDM. In this paper, we consider two interacting forms: 
Q = (3Hpc (denoted as Qi) and Q = l3Hp^^ (denoted as Q 2 )-, where /3 is the dimensionless 
coupling. For convenience, the lACDM models with Qi and Q 2 are denoted as the lACDMl 
and IACDM2 models, respectively, and the ItcCDM models with Qi and Q 2 are denoted as 
the IrcCDMl and Ir(;CDM2 models, respectively. 

This paper is organized as follows. In Sec. 2, we briefly review the two diagnostics, 
i.e., the statefinder hierarchy and the growth rate of structure. In Sec. 3, we explore the 
interacting dark energy models with the geometrical and structural diagnostics. Conclusion 
is given in Sec. 4. 

2 Statefinder hierarchy and growth rate of strnctnre 

In this section, we will first review the statefinder hierarchy diagnostic and then describe 
the growth rate of structure in the interacting DE models. In particular, we present the 
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“parametrized post-Friedmann” framework for interacting DE models, which is used to cal¬ 
culate the growth rate of structure in the lACDM and IrcCDM models. 


2.1 The statefinder hierarchy 

The scale factor of the universe, a{t)/aQ = (1 -|- z)~^, can be Taylor expanded around the 
present epoch to as follows: 

^ = ( 2 . 1 ) 

^ n=l 


where 


Ar. 




ne N, 


( 2 . 2 ) 


with = dAa{t)/. An (n > 3) is transformed to statefinder hierarchy, sli \ in order 

to obtain null diagnostic for the ACDM model, which is expressed as [48]: 


4'^ = ^3, (2.3) 

^ = A 4 3(1 -|- q) (2-4:) 

5 * 5 ^^ = A 5 - 2(4 -I- 3g)(l -I- O'), etc., (2.5) 


where the superscript (1) just is a marked symbol, and q is the deceleration parameter, 

q = -A2. 

In a spatially flat FRW universe consisting of only DE and non-relativistic matter, 
the statefinder hierarchy for ACDM is fixed to be 1 during the cosmic expansion, namely, 
‘S'I^^Iacdm = 1- By using this diagnostic, one can distinguish easily the ACDM model 
from other DE models. However, in this paper, we consider a spatially flat FRW universe 
containing DE, CDM, baryons, and radiation. Obviously, in this case, S'I^^Iacdm is no longer 
a constant, but slowly changes with time during the cosmic expansion. Even so, the ACDM 
model is still treated as a reference model. 

For the interacting DE scenario, we derive the expressions of and in terms of 
the model parameters: 

Sg ^ = 1 -|- -Dde'ic(l w) — —+ 2 D 1 . -|- 2 }{p ’ (^-6) 

Q Q 

^ 4 ^^ =1 — -r(;D^g[3rc(l -|- to) — w'] — -Dde['U’(21 -|- 39t(; -|- ISrc^) -|- Dr(7rc -|- 3rc^ — w') 

- (13 + + 2u,"l - fi,(9 + !2,) - ^(2 + 3w) + ^ + ^. 

In p In p np 

where the prime denotes the derivative with respect to x = Ina and p = The 

above expressions are applicable for all the interacting DE models with arbitrary w{z) and 
Q. For the ItcCDM models, one should let w = const in the above formulae. For the lACDM 
models, one substitutes to = —1 in the above formulae. 
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2.2 The growth rate of structure 

The fractional growth parameter e{z) [50, 51] is defined as 

/(^) 


e(z) = 


/acdm(^) ’ 


( 2 . 8 ) 


where f{z) = dln6/dlna is the growth rate of structure. Here, 6 = dpm/Pm, with p^ and 
6 Pm being the energy density and the density perturbation of matter (including CDM and 
baryons), respectively. 

If the matter density perturbation is linear and there is no interaction between DE and 
CDM, the perturbation equation at late times is written as: 


d + 2H5 = ATrGpmS, 


(2.9) 


where G is the Newton’s gravitational constant. Consequently, the growth rate of the linear 
density perturbation is approximatively given by [55]: 


f{z) ~ ^mizY 


7(z) = 


+ 


3(1 -u;)(l - Iw) 


5-^ 125(1- 


Jl 


(1 - Dm(2)), 


( 2 . 10 ) 

( 2 . 11 ) 


where Djn(^) = pm{z)/3MpH{zY is the fractional density of matter, w is either a constant 
or varies slowly with time. For the ACDM model, 7 ~ 0.55 and e(z) = 1 [55, 56]. For 
other models, the values of e(z) exhibit differences from ACDM, which is the reason why the 
fractional growth parameter e(z) is used as a diagnostic. 

However, for an interacting DE model, the growth rate f{z) cannot be simply param¬ 
eterized by Eqs. (2.10) and (2.11) [57]. For more details of calculating the growth rate in 
interacting DE models, we refer the reader to [57] . But in this paper we calculate the growth 
rate within the “parametrized post-Friedmann” (PPF) framework for interacting dark energy 
[58, 59]. One of the advantages of the PPF approach is that it could resolve the problem 
of early-time super-horizon-scale perturbation divergence in the interacting DE models. We 
will briefly present the PPF framework for interacting DE in the following. 

In an interacting DE model, the conservation laws become 


v.rr = Q^ 


( 2 . 12 ) 


where denotes the energy-momentum transfer of I fluid. Generally, can be split as 


Qj, = a{- Qi{l + AY) - 5QiY, [fi + Qi{v - B)]Y), 


(2.13) 


where 6 Qj and // are the energy transfer perturbation and momentum transfer potential of 
I fluid, respectively, A and B are functions describing the perturbed metric, v denotes the 
velocity perturbation of total matters, and Y and Yi are the eigenfunctions of the Laplace 
operator and its covariant derivatives, respectively. Then, Eqs. (2.12) and (2.13) give the 
following two conservation equations, 


6 p'i + 3{6pj + 6pj) -|- {pj + pi){k}{Vi + 3 H'i) — — AQj), 


(2.14) 


(2.15) 






where kn = k/{Ha) with k the wave number, the prime denotes the derivative with respect to 
X = \B.a, Hl also denotes metric perturbation, and 5pi and 11/ are pressure perturbation and 
anisotropic stress of I fluid, respectively. The values of 5Qi and /i depend on the covariant 
form of Qj. In our work, we choose Qc = —= Q'^c with Uc = a~^ (l — AY, UcT*), so that 
there is no momentum transfer in the CDM frame. Comparing with Eq. (2.13), we can obtain 
^Qde = -5Qc = -PHpc5c for Q = l3Hpc, 6Qde = SQc = -^Hp^eSde for Q = jSHpde, and 
/de = -/c = l3Hpc{v - Uc). 

In above equations. Hi generally vanish for CDM and DE, and the metric perturbations 
are decided by a specific gauge and the Einstein equations. 


Ht B H'j, _ AirGa^ 

^ 3 fc// k'jj k"^ 


<5/9 + 3(p + p) 


V — B 

kn 


A-H' -JL_ 

^ 3 {aHf 


(- 

ykn 


H'j. 

"P” 


AttG V — B 

-{P + P)- 




ku 


(2.16) 

(2.17) 


To complete the perturbation systems of CDM and DE, we also need an extra condition on 
5pi. In uncoupled DE case, one often treats DE as a nonadiabatic fluid and calculates 5pde 
in terms of its rest-frame sound speed, which, however, induces the large-scale instability for 
the interacting DE scenario [60, 61]. To avoid this large-scale instability, we here handle the 
perturbation evolutions using the PPE framework for interacting DE model established in 
[58, 59]. 

Instead of calculating 5pde in terms of rest-frame sound speed of DE, the PPE approach 
completes the perturbation system by establishing a direct relationship between V^e — Vr 
and Vt on the large scales, where Vde and Vt are the velocity perturbation of DE and total 
matters except DE in the comoving gauge. This relationship can be parametrized by a 
function /^(a) as [62, 63] 

lim ^^(pde +Pde)-^%7-^ = -^/c(a)fciffoT- (2.18) 

kH<i H^ kno 


With the help of this extra condition in combination with Eqs. (2.14)-(2.17), we can fi¬ 
nally solve all the perturbation equations and numerically get the growth rate f{z) in the 
interacting DE models. Eor more details, see [58, 59]. 

Besides e{z), we can use a composite null diagnostic (CND), combining e{z) 

with statefinder hierarchy [48]. 

Before we apply the CND to interacting DE models, we wish to make some discussions 
on the consistency relation between expansion history and growth of structure for interact¬ 
ing DE. As is well known, for the non-interacting DE, the matter perturbations are entirely 
determined by the expansion background, as described by Eq. (2.9). This provides a general 
strategy to distinguish between DE and modified gravity (MG) scenarios: a mismatch (dis¬ 
tinct tension) in DE parameter space for constraints from geometrical measurements (e.g., 
luminosity distance and angular diameter distance) and growth-of-structure measurements 
(e.g., linear growth rate) points to a modification of GR. In general, for MG, the growth 
of structure is usually scale dependent, which is a feature not seen in smooth DE models.^ 

^For the MG scenario, in the quasi-static Newtonian regime, we have S + 2H5 — = 0. From 

this equation one sees how the combination of G and >l>/4' alters the linear matter perturbation. If these 
parameters have a scale dependence, then the linear growth function becomes scale dependent. 
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Now, a question arises: Is the growth of structure in the interacting DE scenario also scale 
dependent? Actually, this issue has been addressed in [57] and the answer is “NO”. 

Here we quote the results of [57]. For the = (3HpcUc model, the second-order 
evolution equation for 6 c is 


6 c + 2H6c — 4 :TtG{pc6c + /0b<5b) — 0. (2-19) 

Though this equation has the same form as that for the non-interacting case, the solutions 
6 c are different because the background terms H and pc evolve differently. For the = 
/3HpdeUc model, the evolution equation of 6 c gets two apparent modifications: a modified 
Hubble friction term and a modified source term, expressed as a modified effective Newton 
constant for CDM: 

+ 2h(^1 + Pc - AttGcsPcSc - 47rGpb<5b = 0, (2.20) 

where 

= + ( 2 . 21 ) 
G 3a Pc I V Pc J} 

We see clearly from Eqs. (2.19)-(2.21) that, for these interacting DE models, the growth of 
structure is still scale-independent. Also, from these equations, we find that though the direct 
non-gravitational interaction between DE and CDM is introduced, the matter perturbations 
could still be derived from background observables, such as luminosity distance and angular 
diameter distance data. 

In [64], Alam, Sahni and Starobinsky proposed to reconstruct cosmological matter per¬ 
turbations using observed luminosity distance and angular diameter distance data. They 
considered non-interacting DE cases and thus used Eq. (2.9) to perform this reconstruction. 
Here we have shown that, if one wishes to perform such a reconstruction for the interact¬ 
ing DE cases (ItcCDMl and Ir(;CDM2), one should use Eqs. (2.19)-(2.21) (note that the 
evolution equation for db is not changed and thus is not given here). If the results of these 
reconstructions for both non-interacting DE and interacting DE cases are in tension with 
the actual observational data of growth rate of structure from, e.g., redshift space distortions 
(RSD), then this would give a hint for the MG scenario. 

So far, we have clarified that in the interacting DE models the growth of structure is still 
consistent with the expansion history. Hence, in these cases, both the statefinder hierarchy 
and fractional growth parameter e could be inferred from the same observational data, e.g., 
supernovae and BAG observations. 

3 Exploration of interacting dark energy with geometrical and structural 
diagnostics 

In this section, we explore the deviations from ACDM in the interacting DE models by 
using the diagnostics of statefinder hierarchy and fractional growth parameter e(z). 
Since we only care about the impacts of the EoS parameter w and the coupling parameter 
/3, we fix the other model parameters. For all models, the present-day fractional density 
parameters of CDM, baryons, and radiation are fixed to be: DcO = 0.23, Dbo = 0.04, and 
DfO = 2.469 X 10“^/i“^(I -|- 0.227IAefr) (with the Hubble constant h = 0.7 and the effective 
number of neutrino species Aeff = 3.046). For properly choosing the values of EoS and 
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coupling, we refer to the literature [65-68] . We will first test the impact of w in the IrcCDM 
models. In this test, we fix the coupling to be /3 = 0.02 and consider the cases of w = —0.9, 
— 1, and —1.1. Then, we wish to test the impact of /3 in the lACDM models. In this test 
(with re = —1), we consider the cases of /3 = —0.04, —0.02, 0, 0.02, and 0.04. 

3.1 Testing the impacts of w and /3 in the ItcCDM models 

In this subsection, we will have a look at the impacts of w and j3 in the IrcCDM models using 
the diagnostics of statefinder hierarchy and fractional growth parameter. 


Q ,= / 3 Up ^ 

\ IwCDMI ;5 = 0.02 - w = -0.9 

- w = -1.1 ■ 

IACDM1 /3 = 0.02 • • • •W = -1 
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Figure 1. Evolutions of versus redshift 2 for the IwCDM model with fixed coupling /3 = 0.02. 
The (z) curve of the ACDM model is also shown for comparison. The left panel is for the model 
with Qi = pH Pc and the right panel is for the model with Q 2 = PHpde- 
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Figure 2. Evolutions of e versus redshift z for the IwCDM model with fixed coupling /3 = 0.02. The 
e(z) curve of the ACDM model is also shown for comparison. The left panel is for the model with 
Qi = PHpc and the right panel is for the model with Q 2 = PHp^e- 

First, let us diagnose the IrcCDM models with Qi and Q 2 in the light of the statefinder 
parameter In Fig. 1, we plot the evolutions of versus redshift z for the two IrcCDM 
models with fixed coupling (3 = 0.02 and varying EoS values of ic = —0.9, —1, and —1.1. 
Note that lACDM is a specific case for IrcCDM with w = —1, as labeled in Fig. 1. Since we 
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EoS 

w 



w = —1 



w = —0.9 

w = —1.1 

coupling 

13 

0.04 

0.02 

0 

-0.02 

-0.04 

0.02 

0.02 


c(l) 

0.986 

0.993 

1.000 

1.007 

1.014 

0.698 

1.354 

Qi 

0 ( 1 ) 

‘^4to 

1.032 

1.016 

0.999 

0.983 

0.966 




eo 

0.983 

0.988 

1 

0.998 

1.003 

0.983 

0.992 


c<(l) 

“^Sto 

0.956 

0.978 

1.000 

1.022 

1.044 

0.685 

1.337 

Q 2 

0 ( 1 ) 

‘^4to 

0.975 

0.987 

0.999 

1.013 

1.027 




eo 

0.782 

0.887 

1 

1.098 

1.204 

0.886 

0.887 


Table 1. The current values of the statefinder hierarchy and the fractional growth parameter, 
and Co, for the interacting DE models. 


wish to explore the deviations from the reference model, ACDM, we also show the ACDM 
model as a specific case with /3 = 0 and tc = — 1 in this figure for a direct comparison. We 
find that using only the can easily distinguish the different cases of IrcCDM with fixed 
j3 and varying w. For both ItcCDMl and It(;CDM2, we find that the departures between 
the three lines (with /3 = 0.02) are prominent in the low redshifts up to about z = 3. But 
when we compare the two lines with the same w (i.e., w = —1) and different (3 (i.e., /? = 0 
and 0.02), we find that they are nearly degenerate in the whole redshift range. That is to 
say, it is difficult to probe the deviation of lACDM from ACDM using only the statefinder 
diagnostic We will use the higher hierarchy of statefinder, to further probe the 

deviation from ACDM in the next subsection. 

Next, we diagnose the IruCDM models in the light of the fractional growth parameter e. 
In Fig. 2, we plot the evolutions of e versus 2 ; for the two ItcCDM models where the parameter 
values are taken to be the same with those in Fig. 1. We find that using only the structure 
growth diagnostic, it is also easy to distinguish the cases of IrcCDM with fixed /3 and varying 
w. We find that for both Qi and Q 2 models the e values in the redshift z = 0 are nearly 
degenerate, but they are distinguished well in higher redshifts, with the best distinguishing 
window of z ~ 0.5-1.5. Furthermore, we compare the lACDM and the ACDM cases, i.e., 
the lines with f3 = 0.02 and 0 in this figure. We find that the difference of them is bigger 
in the lower redshifts, and the IACDM2 case (right panel) is better than the lACDMl case 
(left panel). 

In both Figs. 1 and 2, we find that there exist degeneracies in the high-redshift region, 
no matter which diagnostic is used. But this does not influence our diagnosis, since the 
observational data are mainly within the low-redshift region. So we pay more attention to 
distinguish models in the low-redshift region. The current values of and e, i.e., and 
eo, can be viewed as discriminators for testing various cosmological models. For the different 
interacting DE models, the values of and eo are listed in Table 1, which can supply 

assistant information. 

From the above analysis, we find that the impact of w is much stronger than that of 
f3 in the ItcCDM models in both geometrical and structure growth diagnostics. Actually, 
we are more interested in the one-parameter extension to ACDM, i.e., the lACDM scenario. 
Thus, in the next subsection, we will probe the deviation of lACDM from ACDM regarding 
the extra parameter /3, by using the combination of geometrical and structural diagnostics. 
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3.2 Probing deviations from ACDM in the lACDM models 

The lACDM is a one-parameter extension to the ACDM, with the extra parameter /3. We are 
indeed interested in probing the deviation from ACDM in terms of the coupling p. So in this 
subsection we will have a closer look at the lACDM models with the statefinder hierarchy 
and the growth rate of structure. 
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Figure 3. Evolution curves of S' 3 ^^(z), and e{z), and evolutionary trajectories of S' 3 ^^(e) for 

the lACDMl model (Qi = PHpc) with /3 = 0 (ACDM), ±0.02, and ±0.04. In the S' 3 ^^-e plane, the 
current values of {S^^\ e} of the lACDMl model are marked by the round dots. The (e) trajectory 
of the ACDM model is approximately taken to be a point {1,1}, shown as a star. The arrows indicate 
the evolution directions of the (e) trajectories. 

In this case, the deviation from ACDM only comes from the impact of the coupling /3. 
To perform a diagnostic analysis, we take /3 = 0 (ACDM), ±0.02, and ±0.04, for the lACDM 
scenario. For the statefinder hierarchy diagnostic, we consider the evolutions of both S^\z) 
and S' 4 ^^( 2 ;). For the structure growth diagnostic, we consider the evolution of e{z). Moreover, 
we also perform a CND {5'^, ej for the models. 

We first discuss the lACDMl model (with Qi = fiHpc). In Fig. 3, we plot the evolution 
curves of ^^^^(z), 54^^(z), and e{z), as well as the evolutionary trajectories in the plane, 

for the model. In the plane, the current values of ej of the lACDMl model are 

marked by the round dots, and the arrows indicate the evolntion directions of various cases 
for the model. It should be noticed that ~ 1 and e = 1 for the ACDM model during 
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Figure 4. Evolution curves of S' 3 ^^( 2 :), S' 4 ^^(z), and e(z), and evolutionary trajectories of S' 3 ^^(e) 
for the IACDM2 model (Q 2 = PHp^e) with (3 = 0 (ACDM), ±0.02, and ±0.04. In the plane, 

the current values of {S'g^^e} of the IACDM2 model are marked by the round dots. The 
trajectory of the ACDM model is approximately taken to be a point {1,1}, shown as a star. The 
arrows indicate the evolution directions of the (e) trajectories. 


the whole evolution history. As a result, the trajectory of ACDM is a tiny line segment 
in the plane. For simplification, we replace the short S 3 ^^(e) trajectory of ACDM 

approximately with a point {1,1} shown as a star. We find that, for the curves of 
the separations between them are rather near, indicating a near degeneration. But for the 
curves of 54 ^^ ( 2 ), it is easy to hnd that the separations between them become much farther 
than in the case. Specifically, the maximum deviations from ACDM for S^\z) and 

5 * 4 ^^( 2 ) are about 5% and 15%, respectively, showing that could well measure the 

deviation from ACDM and differentiate the geometrical degeneracies in the model. In the 
e(z) panel, we find that the maximum deviation from ACDM for e{z) is only about 2%. Thus 
the e diagnostic only provides tiny differentiation. In the plane, the departures from 

ACDM, from the points of view of both expansion history and growth of structure, could be 
directly measured in this plane. The trajectories also exhibit distinctive features for 

the various cases of the lACDMl model. 

Next, we discuss the IACDM2 model (with Q 2 = (3Hpc). The evolution curves of 53 ^^(z), 
5 * 4 ^^ (z), and e{z), and the evolutionary trajectories in the plane for the IACDM2 model. 
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are plotted in Fig. 4. When we look at the curves, we hnd that they are highly 

degenerate in the high redshift region, but they can separate from each other in the low 
redshift region. Even though in the low redshift range the 53 ^^ (z) curves can be discriminated 
to some extent, the maximum deviations from ACDM for 53 ^^(z) and are only about 

4% and 3%, respectively. Comparing to the I ACDM 1 case (notice the different scales in Figs. 
3 and 4 and refer also to the concrete values in Table 1), we find that actually this statefinder 
diagnostic only provides a mild differentiation and the curves are degenerate in some degree. 
Interestingly, we find that in this model the (S' 4 ^^( 2 ;) curves exhibit even stronger degeneracy 
than the case of as shown in Fig. 4. Thus, for the IACDM2 model, the 5 * 4 ^^ diagnostic 

does not provide any useful help for differentiating the geometrical degeneracies. However, 
we find that for this model the e(z) diagnostic could well discriminate the various cases of the 
model and finely measure the deviations from in the low redshift region (the maximum 
deviation from ACDM for e{z) is about 20%). In the high redshift region, the IACDM2 
model is highly degenerate with ACDM. Therefore, for the IACDM2 model, combining the 
diagnostics of expansion history and growth of structure becomes extremely important. We 
thus perform a CND e} for the model. We find that in the plane the trajectories 

exhibit distinctive features and the departures from ACDM could be well measured. The 
power of CND is clearly shown in this case. 

In the above discussions, the deviation from ACDM for each diagnostic parameter has 
been measured, theoretically. Now we wish to compare these deviations with the uncertainties 
on diagnostic parameters given by observations. To do this, we should first estimate the 
parameter space of the diagnostic parameters using observational data. However, this kind 
of work needs a series of calculations that are not the main focus of this paper. So here 
we quote the observed uncertainties on diagnostic parameters obtained in previous works. 
Uncertainties on and 5 * 4 ^^ from current observations or estimated from simulations can be 
found in [34, 69-74] . It was found that they are rather difficult to measure—for most of the fit 
results, the uncertainties on and at z = 0 are not less than 6 % and 20 %, respectively. 
These uncertainties on and 54 ^^ are larger than the maximum differences [5% (4%) and 
15% (3%)] between the ACDM model and the lACDMl (IACDM2) model, indicating that 
neither the lACDMl model nor the IACDM2 model can be well distinguished from the 
ACDM model. On the other hand, the fit results of e can be found in [50, 51, 64], where the 
uncertainty on e can be constrained to about 5%. With this precision, the IACDM2 model, 
where the deviation from the ACDM model for e is about 20%, can be well distinguished. 
However, the lACDMl model still cannot be distinguished (the deviation for e in this model 
is about 2%). In summary, the diagnostic parameters and e still do not have sufficient 
precisions to completely distinguish lACDM models from the ACDM model, currently. Note 
that their precisions may be further degraded, if one simultaneously constrains S'n^^ and e 
using the same observations, due to the correlations between and e. This situation 
somehow weakens the power of CND {S'g^^e} considered in this paper. Nevertheless, the 
discussion of this paper is still meaningful, as it points out what precisions of the diagnostic 
parameters should be achieved to distinguish the lACDM models from the ACDM model. We 
expect that future observations could measure the diagnostic parameters with such precisions. 
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4 Conclusion 


Currently, the ACDM cosmology provides an excellent description of various cosmologi- 
cal/astrophysical observations. Thus, any new physics beyond the ACDM cosmology, if 
detected and confirmed, would be viewed as a major breakthrough in cosmology and fun¬ 
damental physics. The simplest extensions to the base ACDM cosmology regarding DE 
component includes rcCDM and lACDM. The former considers the model in which the cos¬ 
mological constant A is replaced by some DE held/fluid with constant re, and the latter 
considers the model in which the vacuum energy (A) directly couples to CDM in some phys¬ 
ically profound way, and both are one-parameter extension to ACDM. The mixture of the 
two is called the IrcCDM model. 

In this paper, our focus is the lACDM cosmology. We wish to explore the deviation of 
lACDM from ACDM regarding the extra coupling parameter (3 by using the diagnostics of 
both statefinder hierarchy and growth rate of structure. But, to make our discussion more 
generic, we begin with the more general model, namely, the ItcCDM scenario. We consider 
two interacting form in this paper: Qi = j3Hpc and <52 = PHpde- For the statefinder 
hierarchy, we derived the analytical expressions of and in terms of cosmological 
parameters for interacting dark energy. For the growth rate of structure, we employed the 
PPF theoretical framework for interacting dark energy to numerically obtain e{z) for the 
models considered in this paper. We tested the IrcCDM models with Qi and Q 2 using 
the and e diagnostics. We found that in both the geometrical and structure growth 
diagnostics the impact of w is much stronger than that of /3 in the ItcCDM models. 

We wish to have a closer look at the lACDM scenario that is a one-parameter extension 
to ACDM with the statefinder hierarchy and the growth rate of structure. We plotted the 
evolution curves of S' 4 ^^( 2 :), and €(z), as well as the evolutionary trajectories in the 

Sg^^-e plane, for the two lACDM models with Qi and Q 2 - We found that, for the model with 
Qi, the geometrical diagnostic 5*4^^ could well measure the deviation from ACDM and finely 
differentiate degeneracies in the model. But for the model with Q 2 , we found that neither 
Sg^^ nor ^ 4 ^^ diagnostic could provide useful help in discriminating various degenerate cases 
in the model, while the e diagnostic is fairly helpful for the model. We found that for both 
models the evolutionary trajectories in the ^Sg^^-e plane exhibit distinctive features and the 
departures from ACDM could be well measured. Thus we have shown that the composite 
null diagnostic {S'g^^e} is a promising tool for diagnosing the lACDM models. Besides, we 
also compare our results with the observed uncertainties on diagnostic parameters. We found 
that the diagnostic parameters still do not have sufficient precisions to completely distinguish 
lACDM models from the ACDM model, currently. Anyway, our work gives what precisions of 
measurements should be achieved to distinguish the lACDM models from the ACDM model. 

Acknowledgments 

This work was supported by the National Natural Science Foundation of China under Grant 
No. 11175042 and the Fundamental Research Funds for the Central Universities under Grant 
No. N140505002 and No. N140506002. 


- 12 - 



References 


[1] A. G. Riess et al. [Supernova Search Team Collaboration], “Observational evidence from 
supernovae for an accelerating universe and a cosmological constant,” Astron. J. 116, 1009 (1998) 
[astro-ph/9805201]. 

[2] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], “Measurements of Omega and 
Lambda from 42 high redshift supernovae,” Astrophys. J. 517, 565 (1999) [astro-ph/9812133]. 

[3] M. Tegmark et al. [SDSS Collaboration], “Cosmological parameters from SDSS and WMAP,” 
Phys. Rev. D 69, 103501 (2004) [astro-ph/0310723]. 

[4] D. N. Spergel et al. [WMAP Collaboration], “First year Wilkinson Microwave Anisotropy Probe 
(WMAP) observations: Determination of cosmological parameters,” Astrophys. J. Suppl. 148, 

175 (2003) [astro-ph/0302209]. 

[5] P. .J .E. Peebles and B. Ratra “Cosmology With A Time Variable Cosmological ’Constant’,” 
Astrophys. J. 325 L17 (1988). 

[6] B. Ratra and P. J. .E Peebles “Cosmological Consequences Of A Rolling Homogeneous Scalar 
Field,” Phys. Rev. D 37 3406 (1988). 

[7] R. R. Caldwell, “A Phantom menace?” Phys. Lett. B 545, 23 (2002) [astro-ph/9908168]. 

[8] X. Zhang, “An interacting two-fluid scenario for quintom dark energy,” Commun. Theor. Phys. 
44, 762 (2005). 

[9] X. Zhang, “Coupled quintessence in a power-law case and the cosmic coincidence problem,” 

Mod. Phys. Lett. A 20, 2575 (2005) [astro-ph/0503072]. 

[10] X. Zhang, “Reconstructing holographic quintessence,” Phys. Lett. B 648, 1 (2007) 

[astro-ph /0604484]. 

[11] A. Sen, “Tachyon matter,” JHEP 0207, 065 (2002) [hep-th/0203265]. 

[12] J. Zhang, X. Zhang and H. Liu, “Holographic tachyon model,” Phys. Lett. B 651, 84 (2007) 
[arXiv:0706.1185 [astro-ph]]. 

[13] N. Arkani-Hamed, H. -C. Cheng, M. A. Luty and S. Mukohyama, “Ghost condensation and a 
consistent infrared modification of gravity,” JHEP 0405, 074 (2004) [hep-th/0312099]. 

[14] X. Zhang, “Dynamical vacuum energy, holographic quintom, and the reconstruction of 
scalar-field dark energy,” Phys. Rev. D 74, 103505 (2006) [astro-ph/0609699]. 

[15] J. F. Zhang, X. Zhang and H. Y. Liu, “Reconstructing generalized ghost condensate model 
with dynamical dark energy parametrizations and observational datasets,” Mod. Phys. Lett. A 
23, 139 (2008) [astro-ph/0612642]. 

[16] M. Li, “A model of holographic dark energy,” Phys. Lett. B 603, 1 (2004) [hep-th/0403127]. 

[17] M. Li and R. -X. Miao, “A New Model of Holographic Dark Energy with Action Principle,” 
arXiv:1210.0966 [hep-th]. 

[18] H. Wei and R. -G. Cai, “A New Model of Agegraphic Dark Energy,” Phys. Lett. B 660, 113 
(2008) [arXiv:0708.0884 [astro-ph]]. 

[19] C. Gao, F. Q. Wu, X. Chen and Y. -G. Shen, “A Holographic Dark Energy Model from Ricci 
Scalar Curvature,” Phys. Rev. D 79, 043511 (2009) [arXiv:0712.1394 [astro-ph]]. 

[20] J. F. Zhang, X. Zhang and H. Y. Liu, “Holographic dark energy in a cyclic universe,” Eur. 
Phys. J. C 52, 693 (2007) [arXiv:0708.3121 [hep-th]]. 

[21] X. Zhang, “Heal the world: Avoiding the cosmic doomsday in the holographic dark energy 
model,” Phys. Lett. B 683, 81 (2010) [arXiv:0909.4940 [gr-qc]]. 


13 - 



[22] J. F. Zhang, M. M. Zhao, J. L. Cui and X. Zhang, “Revisiting the holographic dark energy in a 
non-flat universe: alternative model and cosmological parameter constraints,” Eur. Phys. J. C 74, 
no. 11, 3178 (2014) [arXiv: 1409.6078 [astro-ph.CO]]. 

[23] A. De Felice and S. Tsujikawa, “f(R) theories,” Living Rev. Rel. 13, 3 (2010) [arXiv:1002.4928 
[gr-qc]]. 

[24] S. Nojiri and S. D. Odintsov, “Dark energy, inflation and dark matter from modified F(R) 
gravity,” TSPU Bulletin N 8(110), 7 (2011) [arXiv:0807.0685 [hep-th]]. 

[25] G. R. Dvali, G. Gabadadze and M. Porrati, “4-D gravity on a brane in 5-D Minkowski space,” 
Phys. Lett. B 485, 208 (2000) [hep-th/0005016]. 

[26] V. Sahni and A. A. Starobinsky, “The Gase for a positive cosmological Lambda term,” Int. J. 
Mod. Phys. D 9, 373 (2000) [astro-ph/9904398]. 

[27] T. Padmanabhan, “Gosmological constant: The Weight of the vacuum,” Phys. Rept. 380, 235 
(2003) [hep-th/0212290]. 

[28] R. Bean, S. M. Garroll and M. Trodden, “Insights into dark energy: interplay between theory 
and observation,” astro-ph/0510059. 

[29] E. J. Gopeland, M. Sami and S. Tsujikawa, “Dynamics of dark energy,” Int. J. Mod. Phys. D 
15, 1753 (2006) [hep-th/0603057]. 

[30] V. Sahni and A. Starobinsky, “Reconstructing Dark Energy,” Int. J. Mod. Phys. D 15, 2105 
(2006) [astro-ph/0610026]. 

[31] M. Kamionkowski, “Dark Matter and Dark Energy,” arXiv:0706.2986 [astro-ph]. 

[32] M. Li, X. D. Li, S. Wang and Y. Wang, “Dark Energy,” Commun. Theor. Phys. 56, 525 (2011) 
[arXiv: 1103.5870 [astro-ph.GO]]. 

[33] V. Sahni, T. D. Saini, A. A. Starobinsky and U. Alam, “Statefinder: A New geometrical 
diagnostic of dark energy,” JETP Lett. 77, 201 (2003) [Pisma Zh. Eksp. Teor. Fiz. 77, 249 
(2003)] [astro-ph/0201498]. 

[34] U. Alam, V. Sahni, T. D. Saini and A. A. Starobinsky, “Exploring the expanding universe and 
dark energy using the Statefinder diagnostic,” Mon. Not. Roy. Astron. Soc. 344, 1057 (2003) 
[astro-ph/0303009]. 

[35] V. Sahni, A. Shafieloo and A. A. Starobinsky, “Two new diagnostics of dark energy,” Phys. 

Rev. D 78, 103502 (2008) [arXiv:0807.3548 [astro-ph]]. 

[36] A. Shafieloo, V. Sahni and A. A. Starobinsky, “A new null diagnostic customized for 
reconstructing the properties of dark energy from BAG data,” Phys. Rev. D 86, 103527 (2012) 
[arXiv:1205.2870 [astro-ph.GO]]. 

[37] V. Sahni, A. Shafieloo and A. A. Starobinsky, “Model independent evidence for dark energy 
evolution from Baryon Acoustic Oscillations,” Astrophys. J. 793, no. 2, L40 (2014) 

[arXiv: 1406.2209 [astro-ph.GO]]. 

[38] X. Zhang, “Statefinder diagnostic for coupled quintessence,” Phys. Lett. B 611, 1 (2005) 
[astro-ph/0503075]. 

[39] X. Zhang, “Statefinder diagnostic for holographic dark energy model,” Int. J. Mod. Phys. D 14, 
1597 (2005) [astro-ph/0504586]. 

[40] X. Zhang, F. Q. Wn and J. Zhang, “New generalized Chaplygin gas as a scheme for unification 
of dark energy and dark matter,” JCAP 0601, 003 (2006) [astro-ph/0411221]. 

[41] J. Zhang, H. Liu and X. Zhang, “Statefinder diagnosis for the interacting model of holographic 
dark energy,” Phys. Lett. B 659, 26 (2008) [arXiv:0705.4145 [astro-ph]]. 


- 14 - 



[42] M. R. Setare, J. Zhang and X. Zhang, “Statefinder diagnosis in non-flat universe and the 
holographic model of dark energy,” JCAP 0703, 007 (2007) [gr-qc/0611084]. 

[43] C. -J. Feng, “Statefinder Diagnosis for Ricci Dark Energy,” Phys. Lett. B 670, 231 (2008) 
[arXiv:0809.2502 [hep-th]]. 

[44] M. -L. Tong and Y. Zhang, “Cosmic age, Statefinder and Om diagnostics in the decaying 
vacuum cosmology,” Phys. Rev. D 80, 023503 (2009) [arXiv:0906.3646 [gr-qc]]. 

[45] L. Zhang, J. Cui, J. Zhang and X. Zhang, “Interacting model of new agegraphic dark energy: 
Cosmological evolution and statefinder diagnostic,” Int. J. Mod. Phys. D 19, 21 (2010) 
[arXiv:0911.2838 [astro-ph.CO]]. 

[46] F. Yu and J. -F. Zhang, “Statefinder diagnosis for the extended holographic Ricci dark energy 
model without and with interaction,” Commun. Theor. Phys. 59, 243 (2013) [arXiv: 1305.2792 
[astro-ph.CO]]. 

[47] J. -L. Cui and J. -F. Zhang, “Comparing holographic dark energy models with statefinder,” 
Eur. Phys. J. C 74, 2849 (2014) [arXiv: 1402.1829 [astro-ph.CO]]. 

[48] M. Arabsalmani and V. Sahni, “The Statefinder hierarchy: An extended null diagnostic for 
concordance cosmology,” Phys. Rev. D 83, 043501 (2011) [arXiv:1101.3436 [astro-ph.CO]]. 

[49] J. E. Zhang, J. L. Cui and X. Zhang, “Diagnosing holographic dark energy models with 
statefinder hierarchy,” Eur. Phys. J. C 74, no. 10, 3100 (2014) [arXiv: 1409.6562 [astro-ph.CO]]. 

[50] V. Acquaviva, A. Hajian, D. N. Spergel and S. Das, “Next Generation Redshift Surveys and 
the Origin of Cosmic Acceleration,” Phys. Rev. D 78, 043514 (2008) [arXiv:0803.2236 [astro-ph]]; 

[51] V. Acquaviva and E. Gawiser, “How to Falsify the GR-|-LambdaCDM Model with Galaxy 
Redshift Surveys,” Phys. Rev. D 82, 082001 (2010) [arXiv: 1008.3392 [astro-ph.CO]]. 

[52] R. Myrzakulov and M. Shahalam, “Statefinder hierarchy of bimetric and galileon models for 
concordance cosmology,” JCAP 1310, 047 (2013) [arXiv:1303.0194 [gr-qc]]. 

[53] J. Li, R. Yang and B. Chen, “Discriminating dark energy models by using the statefinder 
hierarchy and the growth rate of matter perturbations,” JCAP 1412, no. 12, 043 (2014) 
[arXiv:1406.7514 [gr-qc]]. 

[54] Y. Hu, M. Li, N. Li and S. Wang, “Impacts of different SNLS3 light-curve fitters on 
cosmological consequences of interacting dark energy models,” arXiv: 1501.06962 [astro-ph.CO]. 

[55] L. M. Wang and P. J. Steinhardt, “Cluster abundance constraints on quintessence models,” 
Astrophys. J. 508, 483 (1998) [astro-ph/9804015]. 

[56] E. V. Linder, “Cosmic growth history and expansion history,” Phys. Rev. D 72, 043529 (2005) 
[astro-ph/0507263]. 

[57] G. Caldera-Cabral, R. Maartens and B. M. Schaefer, “The Growth of Structure in Interacting 
Dark Energy Models,” JCAP 0907, 027 (2009) [arXiv:0905.0492 [astro-ph.CO]]. 

[58] Y. H. Li, J. E. Zhang and X. Zhang, “Parametrized Post-Friedmann Framework for Interacting 
Dark Energy,” Phys. Rev. D 90, no. 6, 063005 (2014) [arXiv: 1404.5220 [astro-ph.CO]]. 

[59] Y. H. Li, J. E. Zhang and X. Zhang, “Exploring the full parameter space for an interacting 
dark energy model with recent observations including redshift-space distortions: Application of 
the parametrized post-Friedmann approach,” Phys. Rev. D 90, no. 12, 123007 (2014) 

[arXiv: 1409.7205 [astro-ph.CO]]. 

[60] J. Valiviita, E. Majerotto and R. Maartens, “Instability in interacting dark energy and dark 
matter fluids,” JCAP 0807, 020 (2008) [arXiv:0804.0232 [astro-ph]]. 

[61] J. -H. He, B. Wang and E. Abdalla, “Stability of the curvature perturbation in dark sectors’ 
mutual interacting models,” Phys. Lett. B 671, 139 (2009) [arXiv:0807.3471 [gr-qc]]. 


- 15 - 



[62] W. Hu, “Parametrized Post-Friedmann Signatures of Acceleration in the CMB,” Phys. Rev. D 
77, 103524 (2008) [arXiv:0801.2433 [astro-ph]]. 

[63] W. Fang, W. Hu and A. Lewis, “Crossing the Phantom Divide with Parameterized 
Post-Friedmann Dark Energy,” Phys. Rev. D 78, 087303 (2008) [arXiv:0808.3125 [astro-ph]]. 

[64] U. Alam, V. Sahni and A. A. Starobinsky, “Reconstructing Cosmological Matter Perturbations 
using Standard Candles and Rulers,” Astrophys. J. 704, 1086 (2009) [arXiv:0812.2846 [astro-ph]]. 

[65] J. J. Ceng, J. F. Zhang and X. Zhang, “Quantifying the impact of future Sandage-Loeb test 
data on dark energy constraints,” JCAP 1407, 006 (2014) [arXiv: 1404.5407 [astro-ph.CO]]. 

[66] J. J. Geng, Y. H. Li, J. F. Zhang and X. Zhang, “Redshift drift exploration for interacting dark 
energy,” arXiv:1501.03874 [astro-ph.CO]. 

[67] Y. H. Li and X. Zhang, “Large-scale stable interacting dark energy model: Cosmological 
perturbations and observational constraints,” Phys. Rev. D 89, no. 8, 083009 (2014) 
[arXiv:1312.6328 [astro-ph.CO]]. 

[68] S. Wang, Y. Z. Wang, J. J. Geng and X. Zhang, “Effects of time-varying (3 in SNLS3 on 
constraining interacting dark energy models,” Eur. Phys. J. C 74, no. 11, 3148 (2014) 

[arXiv: 1406.0072 [astro-ph.CO]]. 

[69] E. E. O. Ishida, “Statefinder Revisited,” Braz. J. Phys. 35, 1172 (2005) [astro-ph/0609614]. 

[70] D. Rapetti, S. W. Allen, M. A. Amin and R. D. Blandford, “A kinematical approach to dark 
energy studies,” Mon. Not. Roy. Astron. Soc. 375, 1510 (2007) [astro-ph/0605683]. 

[71] F. Y. Wang, Z. G. Dai and S. Qi, “Probing the cosmographic parameters to distinguish 
between dark energy and modified gravity models,” Astron. Astrophys. 507, 53 (2009) 
[arXiv:0912.5141 [astro-ph.CO]]. 

[72] L. Xu and Y. Wang, “Cosmography: Supernovae Union2, Baryon Acoustic Oscillation, 
Observational Hubble Data and Gamma Ray Bursts,” Phys. Lett. B 702, 114 (2011) 
[arXiv:1009.0963 [astro-ph.CO]]. 

[73] S. Capozziello, R. Lazkoz and V. Salzano, “Comprehensive cosmographic analysis by Markov 
Chain Method,” Phys. Rev. D 84, 124061 (2011) [arXiv: 1104.3096 [astro-ph.CO]]. 

[74] S. Capozziello, M. De Laurentis, O. Luongo and A. Ruggeri, “Cosmographic Constraints and 
Cosmic Fluids,” Galaxies 1, 216 (2013) [arXiv:1312.1825 [gr-qc]]. 


- 16 - 



